Energy-Charge Dependence for Q-Balls 
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o : 

, We show that many numerically established properties of Q-balls can be understood in terms of analytic 
' approximations for a certain type of potential. In particular, we derive an explicit formula between the 
[ energy and the charge of the Q-ball valid for a wide range of the charge Q . 

^ ■ 1 Introduction 
T— I ! 

' As shown by Coleman [1], the existence of Q-baUs is a general feature of scalar field the- 
^ ■ ories carrying a conserved ;7(l) charge [2]. Q-balls can be understood as bound states of 
O ■ scalar particles and appear as stable classical solutions (nontopological solitons) carrying 
^ I a rotating time dependent internal phase. They are characterized by a conserved nontopo- 
^ ' logical charge Q (Noether charge) which is responsible for their stability (see, for example, 
O ■ Refs. [3]- [4]). These features differentiate the Q-ball interaction properties from those 
^ ! of the topological solitons since here the charge Q can take arbitrary values in a specific 
I I range, allowing for the possibility of charge transfer between solitons during the interaction 
^1 process. 

^ ■ The concepts associated with Q-balls are extremely general and occur in a wide 

>■ ! variety of physical contexts [5]. Q-balls are a generic consequence of the Minimal Super- 
! symmetric Standard Model (MSSM) [6] where leptonic and baryonic balls may exist. In 
;_i ' this context the conserved charge is associated with the U{1) symmetries leading to baryon 
. i and lepton number conservation, and the relevant U{1) fields correspond to either squark 
or slepton particles. Thus, the Q-balls can be thought of as condensates of either squark or 
slepton particles. It has been suggested that such condensates can affect baryogenesis via 
the Affleck-Dine mechanism [7] during the post-inffationary period of the early universe. 
Then, two interesting possibilities occur: (i) If the Q-balls are stable and avoid evaporation 
into lighter stable particles like protons, they are cosmologically important since they can 
contribute to the dark matter content of the universe [8]; (ii) If they are unstable, they 
decay in a nontrivial way into baryons protecting them from erasure through sphaleron 
transitions [9]. 

Up till now, comprehensive studies of these objects have been made by using either 
numerical simulations [3]-[ll] or some analytic considerations [1, 12, 13]. In our approach 
we will semi-analytically identify the explicit relation between the energy and the charge 
of the Q-balls by using a semi-Bogomolny argument in the energy density. Subsequently, 
we will show that similar results can be obtained by using the Woods-Saxon ansatz for 
describing the Q-ball profile function; an approach which has been successfully applied 
to describe analytically the properties of multi-skyrmions in three [14] and two spatial 
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dimensions [15]. This way, some universal properties of Q-balls in the thin- wall limit can 
be established. 



2 Energy of the Q-balls 

Although Q-balls can exist in a variety of field theoretical models, we will consider the 
U{1) Goldstone model describing a single complex scalar field <j) in three spatial dimensions 
with potential U{\(j)\). The Lagrangian is 



c = -d^^d'^^-um) (1) 



where the potential is only a function of |0| and has a single minimum at = 0. This is 
equivalent of stating that there is a sector scalar particles (mesons) carrying U{1) charge 
and having mass squared equal to \U"{0). The corresponding energy functional is given by 

E = l (^|5t<A|^ + + UM)^ d^x. (2) 

The model has a global U{1) symmetry leading to the conserved Noether current 

The conserved Noether charge Q is 

Q = ^^J{^dtci)-^dt^)d^x. (4) 

A stationary Q-ball solution has the form 

<t> = e'^'fir) (5) 
where /(r) is a real radial profile function which satisfies the ordinary diff'erential equation 

with boundary conditions /(oo) = and /'(O) = 0. 

This equation can either be interpreted as describing the motion of a point particle 
moving in a potential with friction [1], or in terms of Euclidean bounce solutions [16]. In 
each case the effective potential being C/e//(/) = - U{f) leads to constraints on the 

potential U{f) and the frequency u> in order for a Q-ball solution to exist. Firstly, the 
effective mass of / must be negative. If we consider a potential [/(/) which is non-negative 
and satisfies U{0) = U'{0) = 0, U"{0) = t^;^ > then one can deduce that uj < u;+. Furthermore, 
the minimum of U{f )/p must be attained at some positive value of /, say < /o < oo and 
existence of the solution requires that lj > u>- where = 2U{fo)/fo. Hence, Q-balls exist 
for all u> in the range a;_ < |w| < u+. 

Thus, the charge Q of a stationary Q-ball solution simplifies to 

Q = cjl[f] 

= Attoj J r^f{r)dr (7) 

where /[/] is the moment of inertia. Numerical and analytical methods have shown that 
when the internal frequency is close to the minimal value w_ , the profile function is almost 



constant, implying that the charge (7) is large (thin- wall approximation). On the other 
hand, when the internal frequency approaches the maximal value w+ the profile function 
falls off very quickly (thick- wall approximation). In the thick- wall approximation the 
behavior of the charge Q depends on the particular form of the potential and the number 
of dimensions [13]. In the case studied here we show that Q ^ oo as w ^ 

In order to derive the minimum of the energy of the Q-ball at fixed charge Q, it is 
convenient to represent the energy in the form 

^^« = ^+4./(4! + f/(/)).^dr (8) 

where the stabilizing role of the Q-ball rotational part is obvious. Note that, without 
rotational energy the Q-ball would collapse since i^^ ^ as /(r) everywhere except 
at the origin r = 0. This is similar to the case of rotating skyrmions when the 4-th order 
Skyrme term is omitted in the Lagrangian. The zero mode (rotational) quantum correction 
to the energy, which is proportional to {I[f]y^ plays a stabilizing role in this case. 

The choice of the potential is not unique, the standard requirement is that the 
function U{f)/ p has a local minimum at some value of / different from zero. Here we will 
consider the following potential 

U{f) = f[l + {l- ff] . (9) 
Note that w+ = 2 and = V2 and so that stable Q-balls exist for V2 < w < 2. 



3 Semi-Bogomolny Argument 

We now proceed to obtain an ansatz for the profile function by applying a semi-Bogomolny 
argument [17] in the energy functional (8-9). Initially, this approach was applied to the 
Skyrme model [18], where it was shown that the lower energy bound is proportional to 
the topological charge. Although the model studied here is not a topological one (ie there 
is no topological charge), the Bogomolny argument can still be applied and leads to an 
upper energy bound. This way, the profile function satisfies an exactly soluble first order 
differential equation and the corresponding energy and charge density can be easily derived. 
The same approach was applied for description of multi-skyrmions properties, as presented 
in [14, 15]. 

The energy (8-9) using the Bogomolny argument can be expressed as: 
= (| + i)Q,4./(i(/',Vf(l-/')^)^=* 



LO 1 

2 w 



g + 47r j (J^+f[l-f)^ Pdr-4n J V2ff{l-f)Pdr 

J V2/'/ (1 - f) P dr. (10) 



The equality is satisfied when the total square term is zero which gives the semi- 
Bogomolny equation 

/' = ~V2/(1 - f) . (11) 
The corresponding profile function has the simple form 

/(r) = , ^ > (12) 

yJl + C^ exp(2V2 r) 



which satisfies the boundary conditions /(O) = + and /(oo) = 0, while /'(O) = 
—y/2Cg/{\ + Cg)^/^. Note that the asymptotic value of (12): / ^ cxp(-\/2r) is in agree- 
ment with the equation of motion (6) for = \/2 (ie for large values of Q). In fact, inside 
the Q-ball /' ~ and / ~ 1, while outside the Q-ball /' = and / = and so (12) describes 
accurately the profile function of the Q-ball outside and inside its region where the last 
term in (10) vanishes. However, (12) does not describe accurately the profile function 
of the Q-ball on its surface (the so-called transition region). Although, in the thin- wall 
approximation the analytical and numerical results (presented in Table 1) converge as Q 
increases since the relative contribution of the surface decreases like at large Q. 

For the specific form of the profile function (12) the charge and the energy can be 
evaluated explicitly. We find that 

Qs ^ -^(y ln(CJ + l[ln(CJp + iz3[-C, 

- {^ + i)Q^^{^i + \MC.)f+^^{l + ^yLi.[-C.]) (13) 

where Lin{z) = '^'"^^'■^^ dy and Lii{y) = ln(l - y) is the polylogarithm function. 

Next, the equation for the energy above has to be minimized with respect to 
Cg while Qg is kept constant. We expect the semi-Bogomolny ansatz to be valid only 
when Cg where the initial "velocity" of the trial profile function /(r) tends to zero. In 
this region, the logarithms will dominate the dilogarithm and the trilogarithm functions 
respectively, since these functions tend to zero like polynomials. Thus, by substituting 
z = - In(C^) for z > one obtains 

One can explicitly solve w in terms of the charge given above in order to obtain that 



which when substituted into the energy gives 



(15) 



-KziiT^ + z^) 12^/2 4 V 6 J 

Now, upon minimizing the energy with respect to z we find that the frequency w and the 
charge are given by 



7r2 + 3^2 ' 

2 

2 , ',^2\ 



144 (7r2 + 3^2) 
Finally, the relation between z and w is 

2 



Ql = (6 + 6z + 7r2 + 3z2). (17) 



(.2_2) ^ + V^-12(-^-^) ^^'^^-^^ • 



It is obvious that in the hmit w ^ \/2 the parameter z goes to infinity in consistency with 
the analytical results (14). Thus, the semi-Bogomolny argument is valid in the thin- wall 
approximation. 

The elimination of the z variable from the equations that define and w can be 
easily performed. By letting e = - 2 one gets 



z = 



l+\ 1+e- 



12 



7rV2T7 (2A + 18e + [24 + £ (6 + tt^s)] + f| 



(19) 



(20) 



TT (2 + g) 
\/2£2 



1 + e- 



12 



^ 7r^/2(l + g) ^ 



4 + e 



(21) 



It is obvious that ^ 00 as a; ^ (which corresponds to £ ^ o) while the upper energy 
limit for any value of Q (or £) can be obtained from (21). 

Next, by eliminating e between and we get the Q-ball energy-charge depen- 
dence since 



E^ 



32/3^1/3 5^2/3 



+ 



27/6 ' 211/632/3 

77^/3 (20 - 5477^ + 277r4^ 



TT (4 + 3^2) ^4/3 (i7_2167r2) ^3 
36V2 ^ 2592 2i/6 31/3 '^^ 



1944 25/6 32/3 



o;^/' + o(g-^). 



(22) 



Note that for large values of Q (ie for Q^/^ > l) the expression (22) gives the upper energy 
bound. In Table 1 we compare the results obtained from the semi-Bogomolny argument by 
means of (22) with the ones obtained by solving the full second order differential equations 
(6) numerically, for different values of Q. The agreement is impressive and appears to 
extend far beyond the expected range of validity of the semi-Bogomolny argument. 



Q 


E 




^nurn 




47.691 


95.5007 


87.9945 


1.99750 


0.1530 


25.703 


51.6331 


49.8406 


1.98997 


0.3020 


18.064 


36.4665 


36.0640 


1.97737 


0.4460 


15.168 


30.7595 


30.7182 


1.95959 


0.5798 


14.065 


28.6076 


28.6585 


1.93649 


0.7018 


14.193 


28.8511 


28.8982 


1.90788 


0.8093 


16.268 


32.7688 


32.7587 


1.87350 


0.9000 


19.256 


38.2651 


38.2417 


1.83303 


0.9723 


23.678 


46.2186 


46.2262 


1.78606 


1.0244 


34.910 


65.9541 


66.0215 


1.73205 


1.0555 


61.603 


111.220 


111.482 


1.67033 


1.0654 


149.263 


253.840 


254.515 


1.60000 


1.0564 


722.656 


1140.119 


1141.95 


1.51987 


1.0345 



Table 1: Comparison of the energy given by (22) obtained from the semi-Bogomolny 
argument {e^) with the ones obtained by numerical simulations (£■„„„). 



Note that although the semi-Bogomolny argument and therefore the corresponding 
energy-charge formula (22) holds only for large values of Q the agreement between the 
numerical and analytical results holds for a wide range of Q. However, although the profile 



2000 



1500 
1000. 

E 

500 



— 1 — 

500 



Q 



1000 



1500 



Figure 1: Energy-charge dependence for Q-balls. The upper and lower branches correspond to the thick- wall 
and thin- wall approximations, respectively. 



(12) describes very well the energy of the Q-ball as a function of the charge Q, its value at 
the origin /(O) is always smaller than one and differs considerably from /(0)„„m of Table 1. 

In Figure 1 we plot the energy obtained from the semi-Bogomolny argument (22) 
and from the numerical simulations (-E„„„) for a wide range of Q. Note the impressive 
agreement between the two results. The peculiar behavior of Figure 1 for small Q can be 
explained from Figure 2 where the energy and charge (obtained numerically) as functions of 
w are presented. It is obvious that in a range of ui two different values of energy and charge 
exist. Finally, in Figure 3 we plot the charge obtained from the semi-Bogomolny argument 
(20) and from the numerical simulations ((5„„„) in terms of the parameter a = \/4 - in 
the allowed range < w < 2. 

Rem,ark: The semi-Bogomolny method can be applied for other potentials of poly- 
nomial type. For example, in [12] the potential is of the form 

U{f) = /2 [1 + (1 - (23) 

and the corresponding semi-Bogomolny equation (analogous to (11)) is /' 
its solution ^ 

has the same asymptotic behavior for large r as in (12). 



= -/(!-/) while 

(24) 




Figure 3: The charge Q obtained numerically and semi-analytically (20) as a function of the parameter 
a = \/4 - a;2. 



4 Woods-Saxon Ansatz 



A more general form of the ansatz (12) is widely used in nuclear physics describing nuclear 
matter distribution inside heavy nuclei, or potential of nucleon-nucleus interaction. This 
is the so-called Woods-Saxon distribution and the corresponding profile function is given 

by 

fir) = =■ (25) 

'1 + exp [a (r-r^)] 



In this case, three arbitrary parameters (instead of one) appear in the parametrization, and 
thus, more degrees of freedom exist. Here, corresponds to the radius of the Q-ball and 1/a 
defines the thickness of the shell of the Q-ball. At the origin, the values of the profile func- 
tion and its derivative are: /(O) = C^^ (l + e-'"^o)"^^^ and /'(O) = -aC^<,e-"''« (2 + 2e-"''«)"^^l 
Recall that, due to boundary conditions, /'(O) needs to be very small (in fact, zero). This 
limit can be obtained when the product ar^ is large that is in the thin-wall approximation 
where w ^ V2 (see below). The field-theoretical motivation for the Woods-Saxon ansatz 
was presented in the previous section since for C^^ = 1 and = expi-ar^) the two profile 
functions given by (25) and (12) coincide. 

For Q large, the energy of the Q-balls given by (8-9) can be approximately evaluated 
using the Woods-Saxon ansatz (25). To do so, integrals of the following type are used: 

= f C^wsf'' _^2^^ 

J (l-|-exp[a(r-rQ)])" 

= ^{C„sf''{b'l°+2bll + l'^) (26) 

where n = 1,2,3 and b = ar^ while the moment of inertia becomes: /[/] = 47r/i. In general, 
by letting w = exp[a(r - r^)] while dr = dw/aw, the m-power of the integral J„ is defined as: 

'—dw 27 

z(H-w)" ^ ' 

where wq = exp(-6). After some algebra it can be easily shown that the integrals J™ are 
related via the following recursive relation 

rm _ rm ^ rm-1 (1" "^o)" /9oN 



For Wo < 1 one gets: 



n^b, I'^^b-i, is^b-l, 

h2 ip. ip. 1 



Ifi lyi lyi 

A'^y, /|-|--2L, ^3-y-3i (29) 



where the parameter L = 2(1 - 1/4 + 1/9 - 1/16 + ...) ~ 1.644 does not contribute in the leading 
order expansion of the parameter Q~^/^ (which we assume to be small - see below). Then 
the energy of the Q-ball (8 -9) can be approximated as 



where a = 1 - 36"^ /? = 1 - 7 = 1 + 6Lb-'^ and V = |7rr^. The derivative term of the 
energy: Ederiv = 5/ f^r^dr is going to be neglected (initially). 

The minimization of with respect to iVC"^^ occurs at 

{iVClX.. = ^2(2-2C;£a + C^^/3) ^^^^ 
which determines r^. The corresponding minimum of the energy is 

E^c,V2Q{2-2aCl^+pCl^). (32) 
Further minimization with respect to C^g gives the following value for the energy 

£^^^/2Q(2-^) (33) 

(^ws)L„ = ^IP- large 6 (of order Q^/^), the parameter (C^s)„„„ can be approximated 
by (Ci^s)„i„ ^ y^l + 3/(26) = 1 + 3/(46) and in this limit the energy becomes: 

^;^^^/2Q(l + A) (34) 

where 

Note that equation (34) implies that E^ \/2Q as 6 ^ 00. 

Next the derivative energy contribution is considered: 



^ anC^^ ( n 2r 2L\ 

i;deriv - + — + ^ • (36) 



a 

Taking into account the highest order terms in one obtains 

Further minimization with respect to h gives the energy-charge dependence of the Q-balls 
up to order 0{Q'^/^): 

^» - ^«+(il)"«^" 

= 72 g + 1.35702 g2/3 (38) 



which occurs when 



3V2Q 



1/2 



12 N^/^ 



Q ~ 1.563 Q^/^ (39) 



Note that the 0{Q'^/^) terms of (38) and (22) coincide, and therefore the two analytic 
methods are in good agreement. In addition, from (35) and since h = ar^ one gets that 
a^.^ = 2V2 which is the value obtained from the semi-Bogomolny argument (12). 



To conclude, let us state that in the thin-wall limit the profile function (25) is given 
by /(O) ~ 1 + 3/(45) (or by /(O) ^ 1 + l/(2&) including higher order corrections). Moreover, in 
this limit the Woods-Saxon distribution (25) is close to the exact values of the Q-ball profile 
function (obtained numerically) presented in Table 1. In fact, for Q = 149.26: /(O) = 1.06 
and /(0)„„^ = 1.0564 while for Q = 722.66: /(O) = 1.0356 and /(0)„_ = 1.0345. 

Remark: For the alternative potential given by (23) the semi-Bogomolny function 
(24) can be generalized in a similar way as (25) and the corresponding energy is similar to 
(38) since 



Note that the two expressions differ only by a 2/3 factor before the term 0{Q^/^). 
5 Discussion and Conclusions 

It has been shown that for specific parametrizations of the scalar field a semi-analytic 
treatment for Q-balls exists leading to transparent and simple results. Two kinds of 
approximations have been considered: one based on a semi-Bogomolny argument which 
gives an exponential-step parametrization for the profile function and the Woods-Saxon 
parametrization (the semi-Bogomolny generalization) motivated also by nuclear physics 
experience. 

The agreement of the results obtained using both approximations with the numerical 
ones follows from the fact that the ansatz for the profile function obtained from semi- 
analytic approaches have the correct asymptotic behaviour as Q and r are large. This was 
not the case in the semi-analytic treatment of the Skyrme model in three or two spatial 
dimensions [14, 15]. Although the energies obtained were accurate within 0.5% compared 
to the exact ones for large values of baryon number, the asymptotic behaviour of the profile 
function was incorrect [14]. 

The thickness of the Q-ball surface (ie transition region) where the profile / decreases 
from / ~ 1 to / = can be estimated by 



For Q large (and using the results of section 4), the thickness is independent of the charge 
since t ~ \/2. Thus, the large Q-balls can be visualized as spherically symmetric balls with 
constant internal energy density 



These balls are surrounded by a surface of constant thickness and constant average energy 
density per unit volume since 




(40) 



Pe,v 



~ 2. 



(42) 



leriv 



Pe,s — 




4 



(43) 



in natural units of the model. Therefore it is energetically favorable for small Q-balls to 
fuse into a bigger one since the surface of a single big Q-ball is smaller than the sum of 
surfaces of several smaller Q-balls, for the same value of Q (or with the same total volume). 



Our approach can be extended in lower (and also higher) spatial dimensions in a 
natural way. In particular, in the one-dimensional case the energy-charge dependence is 



E{Q) 




(44) 



which is in a close agreement (within 1%) with the numerical results obtained in [19]. In 
addition, the value of the profile function at the origin is (approximately) given by 



where terms of the form expi-ar^) have been neglected since Q is large (ie »)■ 
In the two-dimensional case the corresponding results are 



The formulas (44), (46) and (38) indicate that the relative contribution of the 

surface energy {Edenv) increases as the dimensionality of space increases. This property of 
Q-balls can be useful in cosmological applications. In some respect Q-balls arc similar to 
the multiskyrmions which correspond to bubbles of matter with universal properties of the 
shell, where the mass and baryon number density is concentrated [14, 15]. 

We conclude by saying that our arguments can be extended to other types of po- 
tentials. Another important issue concerns the semi-analytic identification of the profile 
function satisfying equation (6) for (9). We expect to report on this problem soon. 
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